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Matrix Model Conjecture for Exact BS Periods and 

Nekrasov Functions 

A.Mironov* A.Morozov^ and Sh.Shakirov* 

Abstract 

We give a concise summary of the impressive recent development unifying a number of different funda- 
mental subjects. The quiver Nekrasov functions (generalized hypergeometric series) form a full basis for all 
conformal blocks of the Virasoro algebra and are sufficient to provide the same for some (special) conformal 
' blocks of W-algebras. They can be described in terms of Seiberg-Witten theory, with the SW differential 

given by the 1-point resolvent in the DV phase of the quiver (discrete or conformal) matrix model (/?- 
£N| ' ensemble), dS = ydz + 0(e 2 ) = £ 2p p^ ^\z), where e and (3 are related to the LNS parameters ei and c^. 

This provides explicit formulas for conformal blocks in terms of analytically continued contour integrals and 
q ■ resolves the old puzzle of the free-field description of generic conformal blocks through the Dotsenko-Fateev 

integrals. Most important, this completes the GKMMM description of SW theory in terms of integrability 
theory with the help of exact BS integrals, and provides an extended manifestation of the basic principle 
■ which states that the effective actions are the tau-functions of integrable hierarchies. 

m i 

1. Introduction. A renewed interest to SW theory [1]-[14] and to Nekrasov functions [15]-[19] is caused by 
the AGT conjecture [20]-[31], which connects the subject to the classical field of conformal field theory [32]-[38]. 
Today it is clear that a new unification of the fundamental importance emerges, bringing together at a principally 
new level the CFT, the theory of loop algebras, SW theory, quantization theory, Baxter equations, DV phase 
of matrix models, loop equations, the theory of hypergeometric functions, symmetric groups, Hurwitz theory, 
Kontsevich models and modern combinatorics. This unification is capable to resolve a number of long-standing 
problems in each of the fields. The goal of this paper is to briefly summarize our knowledge about this emerging 
pattern, which is scattered and expressed in length in a number of fresh [20]-[31] and older [39]-[51] papers. The 
main emphasize will be put on the description of Nekrasov functions in terms of SW theory, where the exact 
prepotential is expressed through the exact Bohr-Sommerfeld integrals, and the integrand is provided by the 
1-point function of conformal matrix model (or quiver /3-ensemble) in the Dijkgraaf-Vafa (DV) phase [52]- [55]. 
This conjecture, explicitly formulated in [25] and further investigated [29, 30], makes the picture complete, 
resolves the remaining uncertainties (about the shape of the second deformation) in [56, 57] and finalizes the 
program [3] to reformulate SW theory of [1, 2] in terms of the BS integrals and underlying integrable systems. 
We do not discuss long formulas, checks and even evidence in favor of all these conjectures: all calculations in 
these fields remain long and tedious, and most statements still need to be checked and proved, however, the 
entire picture is starting to get relatively clear. 
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2. Nekrasov functions from BS/SW periods. The central object of emerging unification is the exact 
SW-Nekrasov prepotential T{a\ei 1 1^) — £1^2 log^Ar e fc) which now has a number of different interpretations: 

(A) •Z'iVefc = sum over chains of Young diagrams = generalized hypergeometric series 

(B) Zftek = sum over partitions with Plancherel like measure = discretized matrix model 

(C) %Nek = B = conformal block, depending on a number of external and internal dimensions 

(D) ZN e k = partition function of conformal (3 — ensemble in DV phase = generalized Dotsenko — Fateev integrals 

(E) Z Nek = solution to consistent set of SW equations 

ai = f dS ei<e2 , 
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with deformed SW differential 



dS^ 2 = (y + 0(e 2 ))dz = J2e 2p $ ll \z) 



(1) 



• ■ZjVefc is the generalized hypergeometric series [24], a sum over chains of Young diagrams (integer partitions), 
introduced by N.Nekrasov [15] as expansion of the LNS multiple contour integrals [12] originally obtained by the 
localization (Duistermaat-Heckman) technique [58] from regularized integrals over the ADHM moduli spaces 
[59] of instantons [60] (for their original relation to SW theory and detailed references see [13]). 

• Discrete sums over integer "eigenvalues" are immediately associated with the Nekrasov expansion of Z^ek 
into Young diagrams. They are examples of the character expansions of r- functions [61] and have deep and far- 
going relation to combinatorics of symmetric groups [18], to Hurwitz theory [62] with the cut-and-join operators 
[63], and finally to the matrix models of Kontsevich type [41]. 

• The fact that Z^ e k is the same as the conformal block is the celebrated AGT conjecture [20]-[31], originally 
motivated by the study of 5-brane compactifications on Riemann surfaces [64]- [69]. The AGT conjecture turned 
the Nekrasov functions into a very serious candidate on the role of the first special function of the string theory 
[47] and suggested that they indeed generalize the hypergeometric series p F q in the right direction. The AGT 
conjecture also attracted an attention to the old problem of conformal blocks description through the Dotsenko- 
Fateev integrals [34], which are the crucial element of the free- field description of arbitrary conformal models 
[35, 36]. 

• The Dotsenko-Fateev integrals [34] , the correlators of free fields with insertion of the " screening charges" , 



\ a i=l-' C i I a<b i,a i<j 

were long ago put into the context of matrix models [42]- [44] and are associated with what is called conformal 
or quiver [45] matrix models. They are basically defined by the r.h.s. of (2). Then, the AGT conjecture 
immediately implies that these models have a direct relation to Nekrasov functions, as anticipated in [54], 
carefully formulated in [25] and further discussed in profound papers [29, 30]. 

• These papers, however, stopped short from formulating the exact shape of the Seiberg-Witten equations 
for the exact e-deformed prepotential. The problem is to find the shape of exact SW differential dS, which is 
very well known for ei = e 2 = 0, more or less understood (though still at conjectural level, see [56, 57]) for 
e 2 = and remained a mystery for generic t\ and e 2 . Matrix models considerations provide the answer: dS is 
the DV differential ydz plus higher genus corrections, with 

y 2 = W'(z) 2 + 4pf^(z) (3) 

for the conformal matrix model potential 

k 

W(z) = J2 b -^\og(z- qa ) (4) 

a 

and /? = b 2 /e 2 . In the simplest case of the SU(2) 4-point conformal block k = 4, 91,2,3,4 = 0, 1, quv, 00. The 
higher genus corrections are described by the theory of Dijkgraaf-Vafa phases in matrix models, [52]- [55], which 
remains underdeveloped, but is still capable to provide some explicit checks, see [29, 30]. 
Relations between various parameters are summarized in the following list: 

£ 2 = -ei£2, (5) 
The central charge of the related conformal model is 



l-f^l-e^-^) 2 , i.e. e = e 1+ e 2 =e^-^)=b- £ j 



(6) 



so that the screening charge is simply b = e\, another screening would have b' = e 2 = e 2 /b, but it docs not 
appear in our considerations. The conformal dimension of the vertex operator e la< ^ is 

A a = = (7) 

eie 2 e 2 
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Figure 1: Potential — W(z) for four masses a = (3, 2, 1, 3) situated at positions q = (0.5, 1.5, 2.5, 3). The extrema 
of W(z) are filled with the eigenvalues with filling fractions N\,N 2 , N 3 . 
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Figure 2: Diagram for the five-point conformal block, with the dimensions of external states corresponding to 
the parameters a and the dimensions of intermediate states corresponding to the parameters a\ — Ni — N 2 , 
a 2 = N 2 - N 3 . 



3. DV phase of the /3-ensemble. Of all these subjects the central one for this paper is explicit formula 
(1) for exact dS, and it involves two advanced notions in matrix model theory: the /3-ensembles and the DV 
phases. The /3-ensembles of which the conformal matrix models and Dotsenko-Fateev integrals are examples, 
are not, strictly speaking, matrix models: when the Van-der-Monde determinant A(z) = Hi<j( z * — z j) 1S ra ised 
to the power 2(3 with (3^1, 1/2, 2 they are rather not matrix but eigenvalue integrals. However, the theory of 
matrix models [46] does not really have so much to do with literally matrix integrals: it is rather a theory of 
r-functions, subjected to additional constraints, like string equations and, as a consequence, the set of Virasoro 
or W^-like constraints [39], also known as loop equations [40]. These equations are not very much affected by 
the change of (3, and all the "matrix-model" multi-resolvents p^ m ^ [47] (or their Fourier/Laplace transforms, 
relevant in some combinatorial contexts [70, 48]) are equally well defined for arbitrary (3 (see [71] for one of 
the first reviews of the subject and [72] for a fresh look). In principle, for (3 ^ 1 the terms with half-integer 
p (associated with non-orientable Riemann surfaces) could arise in the large N genus expansion for arbitrary 
potentials [25, 73], but this does not necessarily happen in expansion (1). 

Another ingredient of our description of exact dS is the Dijkgraaf-Vafa (DV) phase: a peculiar phase of 
arbitrary matrix model, describing a multi-cut solution to the Virasoro like constraints where the cuts are, as 
usual, associated with extrema of the action, and the numbers Ni of eigenvalues (filling numbers), ascribed 
to each extremum are kept fixed. This phase is well defined for all values of Ni, not obligatory infinite [47], 
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and it is characterized by the peculiar DV differential ydz, which is essentially the 1-point resolvent at genus 
zero p^ x \z). By itself, it provides the non-deformed, original SW differential. Switching on the t\ = —€2 
deformation modifies it by higher genus corrections, while the orthogonal deformation to e = ei + 62 7^ 
corresponds to making j3 ^ 1. 

The need for the DV phase appears because the number of parameters in the Dotsenko-Fateev integral (2) 
docs not match the number of parameters in the conformal blocks. The conformal block depends on dimensions 
of external legs: these are labeled by parameters a a in (2), on the central charge: this is hidden in b (or (3 
or £i/e 2 ), and on dimensions of intermediate states: there is nothing obvious to encode them in (2). This 
mismatch is the long-standing puzzle in conformal field theory. The way out, implied in (1) and suggested, as 
we understand, in [25] is that the DV phase is additionally labeled by parameters iVj (filling numbers), and 
these are the ones that are missed hidden parameters in (2) and they can be used to account for intermediate 
dimensions in the conformal block. The number of the extra parameters is the number of extrema of the 
potential W(z) in (4), and it is exactly the number of intermediate states in the conformal block. Moreover, 
the intermediate a-parameters are linear combinations of the filling factors: 

ai= I dS €u€2 ~ N I+1 - Nj (8) 

J A, 

There are powerful techniques in matrix model theory to deal with DV phases, based on Givental decom- 
position [50] on the formalism of check-operators [49], and on Whitham dynamics [53]. Here we do not go into 
details of these approaches and just sketch two of the possible lines of reasoning, which seem most promising. 



4. Direct description of DV phase: f3 = 1. As a first illustration, we consider a particular version 
of the Givental decomposition method [50, 54]. In this particular form it is applicable only for (3 = 1, but it 
highlights essential properties of the DV phase. For (3 = 1 the conformal matrix model (2) possesses the simple 
determinant representation 

II j dziHizi-qa)"" Hizi-Zj) 2 =detffy (9) 

i=l^ C ' i,a i<j % ' 3 

which admits a straightforward analytical continuation to large differences N a — Nb without any reference to 
the spectral curve formalism. Here Hij is the so-called matrix of moments: 



fltf(C) = J dz z^Yliz-qa)"' 



(10) 



In the DV phase there is a finite set of allowed contours, in one-to-one correspondence with extrema of the 
potential. One substitutes into (9) a linear combination 

fly=5> 7 fly(C 7 ) (11) 
I 

and picks up the coefficient in front of Y[i u? 1 in dot H , this defines the DV phase partition function. For 
example, in the particular case of the four-point function the Hij{C) is the Eulcr hypergeometric integral 

H t] {C)= J dz z t+ i +ai (z- l) a2 (z-q) a3 (12) 
c 

which is easily expressed through hypergeometric functions: depending on the contour, is equal to either 
1 

/ dz ^+°i(*-l) a3 (*-g) as = ( ~ 1)a2+ "V if ~ Q3 ' -«i-«2-a3-l-i-j \ (13) 
J 1 + ?; + j V -a 3 -a!-i-j H ) K 1 



or 



q 

dz z l+ ^{z - l) a i(z - q)^ - 1 ! Y+i + j { 



(.Ijoa+asgl+i+j+as+a, ^ ^ + • + ■ + j 

2 + a 3 + oti + i + j 



(14) 
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The DV phase partition function is given by 



1 



j> j> -J^ det lujdz z l+1+ai (z - l)" 2 (z - q)° 3 +vjdz z l+ 3 +ai (z - l) tt2 (z - 



q)<*> (15) 



Note that the indices i and j run from to N\ + N 2 — 1. This expression allows a straightforward analytical 
continuation to large Ni — N 2 in each order in quv- The explicit calculation for arbitrary ai,a 2 ,a 3 is rather 
lengthy and will be presented elsewhere. In this paper we consider only the case of a\ = a 2 = a 3 = 0, to 
provide a simple and clear evidence of relation between integral (15) and the 4-point conformal block. 
In the case of vanishing c^, the partition function takes the form 

/du f dv ( u v q i+ J +1 \ 

^f^^{i + J+l + i + J+l) (16) 
From dimension counting, the partition function behaves as 

Zn u nM = const(JVi, N 2 ) ■ q N * ■ (C (N U N 2 ) + Ci(JVi, JV 2 )g + C 2 (N 1 ,N 2 )q 2 + . . . ) (17) 

Thus, the quantities of interest are the normalization constant 

const(JVi, JV 2 ) = F Nl (2N 2 ) F N . 2 (0), 

p 1 (a + N\ N ^ r(£ + l/2)r(L-l/2) / (a + L)(a + 2N L) \ L 
FN{a) -8^{-^r) ll 16 r(L)r(L + l) I L(2N-L) ) ^ 

and the expansion coefficients 

C (N U N 2 ) = 1, 

-"Lb*™, 

1 - 3Ni + 8N$ < ,2 1 ?at2 



(19) 



Note that these coefBcients are not symmetric w.r.t. Ni and N 2 . Actually, there are no reasons to expect such 
a symmetry: the integration contours of the first and of the second type arc different. If one would integrate 
from to qi with i = 1,2, then the answer would be symmetric with respect to changing (iVi,gi) <-> (N 2 ,q 2 ). 
This symmetry gets broken, when one chooses q\ = 1, q 2 = q. 

One observes that the poles can emerge in the iV-paramcteres, in fact they are closely related to the t'Hooft- 
De Wit spurious poles in the theory of unitary matrix integrals [75, 61]. These coefficients need to be compared 
with the first two coefficients of the four-point conformal block [22] 

B< >(Ai,A 2 ,A 3 ,A4,A) = l 
^(Ar, A 2 , A 3 , A 4 , A) = (A + Ar - A 2 )(A + A 3 - A 4 ) 

flgW A 2 , A 3 , A 4 , A) ^ A + Al - A ^ A + ^ - A ^% A3 = A4KA + ^ " ^ + ^ + W 
| [(A 2 + Ai)(2A + 1) + A(A - 1) - 3(A 2 - AQ 2 ] [(A 3 + A 4 )(2A + 1) + A(A - 1) - 3(A 3 - A 4 ) 2 ] 

2(2A + 1) (2A(8A - 5) + (2A + l)c) 

On the matrix model side, putting N\ — n 1 N 2 = —n, one finds 

C (n,-n) = 1, 
d(n, -n) = —, 
„. . n 2 (l-6n 2 + 9n 4 + 8n 6 ) 

Ca( "' ^ = 4(1 - 4n 2 ) 2 (21) 
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On the conformal block side this corresponds to putting Ai = A2 = A3 = L 

£ (0) (0,0,0,0,A) = 1, 

4 1) (0,0,0,0,A) = ^ 

A(1-6A + 9A 2 + 8A 3 ) 
4(1 -4A 2 ) 2 



0, when 



£^(0,0, 0,0, A) 



One can see that 



Ci(n,-n) =Bf (0,0, 0,0, n 2 ), i = 0,1,2,... 



(22) 



(23) 



i.e. the perturbative expansion of the matrix model indeed reproduces correctly the perturbative expansion of 
the conformal block. 

Moreover, this observation can be generalized: on the matrix model side, putting Ni = n + 5, N 2 = —n, one 
finds 

C (n + S, -n) = 1, 
n 2 -5 2 



C 2 (n + 5, -n) 



Ci(n + 5, -n) 

(n 2 - S 2 )(-6n 2 -S 2 + 9n 4 + 3n 2 S 2 + 8n 6 - 8n 4 <5 2 + 1) 



4(1 -4n 2 ) 2 

On the conformal block side, one puts Ai = A 2 — A3 = and gets 

B(°)(0,0,0,A 4 ,A) = 1, 

B«(0,0,0,A 4 ,A) = ^^ 

(A - A 4 )(-8A 2 A 4 - A 4 + 3A 4 A + 8A 3 + 9A 2 + 1 - 6A) 
4(1 - 4A 2 ) 2 

Again 



(24) 



4 2) (0,0,0,A 4 ,A) 



C t {n + 6, -n) = B < £ ) (0,0, 0,5 2 ,n 2 ), 



0,1,2,... 



(25) 



(26) 



This two-parametric equality provides more evidence of the equivalence between conformal block and matrix 
integral. Note that in this example Ni — N 2 — 2n + 5 is not proportional to a = n, the additional constant 5/2 
being removable by a proper choice of the contours. Two possible and natural generalizations arc in order: to 
non- vanishing a\, a 2 and 0:3 (they should match the missing parameters Ai, A 2 and A 3 on the conformal block 
side) and to /3 ^ 1. We now proceed to the second one. 

5. Direct description of DV phase: (3^1. For (3 ^ 1 and vanishing masses the partition function of 
the conformal matrix model takes the form 

N ± 1 N 2 1 N t N 2 

z Ni ,n 3 w, 1) = n / du * n / ^ n>* - u ^ n> - v ^ fi n - v ^ ^ 

»=1 g »=1 i<J i<rj i=lj=l 

No determinant representations are available for this partition function; however, the integrand is just a poly- 
nomial and can be integrated explicitly for each particular N±, N 2 and (3. From dimension counting 

Z nInM = constW(JVi, N 2 ) ■ q PNi(N 2 -i)+N 2 . (cW {NuN2) + ^{N^N^q + C^iN^N^q 2 + . . . ) (28) 

Again, the quantities of interest are the normalization constant and the expansion coefficients. After some 
calculation, one finds the normalization constant 



const^A^, N 2 ) = F^(2(3N 2 ) fg(0) 



(29) 



G 



where 



(0L + l) a (20N-0 + 2-0L) ( 



-0) a (20N-20 + 2-0L) c 
and (a)k = T(a + k)/T(a). We are yet unable to find explicitly; only for = 1 we know that 

MD 1 TT 1fi L r(£ + l/2)r(L-l/2) 

r(L)r(L + i) 

Also one finds the coefficients 

C (JVi,JV 2 ) = 1 



(30) 



(31) 



(32) 



Ci(JVi,JV 2 ) = -0 



N? + 2N 1 N2 + f3-l 



2 2 

Comparison with the corresponding coefficient in the conformal block gives 



0-1 



Ci(n,-n) = 0y + ,., 



-n 



(0,0, 0,0, A) 



This is in accordance with the formula 



a(ei + e 2 - a) 



£ie2 



Indeed, on dimensional grounds a = yln and one obtains 



An(ei + £2 — An) 



£l£2 



0-1 



A system of two equations 



2e l£2 2 ' 



A 2 



2eie 2 



I 
2 



(33) 



(34) 



(35) 



(36) 



(37) 



is satisfied by A = — ei, /? = —e 1 /e 2 (or A = — e 2 , /3 = — e 2 /ei). Thus, we derived "from the first principles" 
the Dotsenko-Fateev integral representation of the conformal block. 

Like in the case of = 1 eq.(34) can be easily extended to include A4 7^ 0: 



Ci (n + 6, -n) = B { 1 } (0, 0, 0, A 4 , A) 



(38) 



Because of our choice of contours we have a = — e\(n + 25), 0:4 = —e\5. 



6. Genus expansion and loop equations. We now turn to a more traditional approach to the genus 
expansion [47], straightforwardly applicable for any values of 0. In the case of DV phases it has certain 
peculiarities (the formalism of check-operators [49] has to be very effective in this case). 

The Ward identities for the conformal matrix model follow as usual [39] from invariance under the shift of 
integration variables Szi ~ with n > — 1. Invariance of the r.h.s. of (2) implies 



(N . 
n t dz * ri( 



N 



Zi-q a T a/£ Il^-^) 2/3 =0, 



(39) 



i<3 



i.e. 



20- 2 E 



n+l _ n+l 



+EE 



Ma 2. 



n+l 



9a 



(40) 
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where fj, a = 2ba a = 2e\a a - Summing up these identities over n with the weights £ ™ 2 , one obtains the Virasoro 
constraints in the form of a loop equation 



In (1) one needs the resolvent [47] 



^ |1) (0 = £ 2 (E^:) = E^? 1) (0 



(41) 



(42) 



p=0 



which satisfies the equation 



p (p h1) (o) 2 + WW 1 *® - no = a - /5) £ 2 ^^ ~ ^V-'%,o 



(43) 



where W"(0 - E a |^ 



AT 



pj |2) (a,6) = £ 2 E 



and 



9a - Zj 



(44) 



(45) 



< ... > c here denotes the connected correlator, and the constants c a can not be defined from the symmetry 
considerations (from the Ward identities). 

Eq.(43) does not define the one-point resolvent directly: the best one can do at once is to find the spherical 
term with p = 0, when the first term in (41) factorizes: 



where 



/ (0) (o = E 



£ - Qa 



1 \ r (( 

Er \ v t- a 



(0) 
a 

£ - <7a 



(46) 



(47) 



After rescaling of p^ ! 1 ) the loop equation (46) can be interpreted in terms of average of the stress tensor 

[40, 43, 76, 20, 29], the term at the r.h.s. corresponds to adding the (^ffi — d 2 cf> in the case of non-unit 

central charge. 

Now solving the quadratic equation (46) , one gets [30] : 



p (0|1) (0 - ~ W ' { ^ +m , v 2 - W(fl a + W (O) (0 

Neglecting the total derivative in (48), one obtains that 

dS(J3 = 1) ~ y(z)dz 



(48) 



(49) 



This SW differential has poles at the points q a and also at £ = oo, and according to the general rules of [2] 
residues at poles are the masses of matter supermultiplets (4 = 2Nc fundamental and k — 3 bifundamental for 
k + 1 > 4) . These masses are linear in a a as a part of the AGT relation [20] . 

In the simplest U(2) case of the 4-point conformal block, i.e. for k + 1 = 4, and for all external ai.2,3.4 = 
(thus W'(z) = 0) and also for e = 



ydz = 



y/ud,Z 



\/z{z- l)(z-q uv ) 



(50) 
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Note that c-parameters in (47) are no obliged to vanish together with fj, a , despite they can seem to do so from 
that formula: the possibility to treat c as completely free parameters and take the limits of this kind is the 
specifics of DV phases. The modulus a ~ y/u <~ N\ — N2. For ydz in the case of W(z) ^ see [30], in the 



inverse limit of pure gauge theory ydz — ^f\J u ~ z — \ becomes the sine-Gordon presymplectic form [3]. In 
the zeroth order in e one obtains the classical contribution to the prcpotcntial: T = to 2 + 0(e), with 



QUV dz 



1 , Jo y/z(z-i)(z- quv ) 6>f (r) 

r= 2^=^—1 — =* quv= em (51) 

J0 y/z(z-l)(z-q uv ) 

which gives the celebrated relation between the bare quv and the instanton-corrected dressed coupling q [77, 
20, 27]. This should be compared with the expression through the Nekrasov functions, 

/ quv a 4 \ a 4 1 
^Nek{e = 0) = lim e!£ 2 log 1 H — ^ + 0(q uv ) ) = quv~r^ 2 + • • • = quva 2 + -Auv^v + ■■■ ( 52 ) 



£2 ^o ^~ b V ' ei e 2 4 a 2 -e 2 ■ -wvjj ^" 4a 2_ e 2 1 ■■• ' 4 

In the case of finite parameter e 2 = 0e 2 + 0(e 4 ), and only the first term here survives in the genus zero 
approximation. Thus, indeed the term with e = in (52) is immediately reproduced from our consideration. 

7. e 2 = and other limits. The interesting limit £2 = with finite e\ was considered in much more detail 
in [57]. This corresponds to keeping e = 0, while ~ e\/e 2 . In this case the r.h.s. in (43) can not be neglected, 
instead one gets a typical Ricatti equation for p^ x \ which should be compared with 

i/j(i)-expf — f x pdx) 

(-efd 2 + e ix + e~ ix - u) V = M ' p 2 +e l p' = V{z) (53) 

of [57]. In the case of the pure gauge U{2) theory V(z) = u — z — -, z = e lx is the sine-Gordon potential, while 
the case Nf ^ is described in terms of XXX (at d = 4), XXZ (at d = 5) and XYZ (at = 6) magnetics 
[8, 66, 74], and for ei 7^ the counterpart of Shroedinger like equations of [57] is a somewhat more complicated 
Baxter equation 

(K(-ieid) + e lxl2 K+(-it x dy xl2 + e -" /2 X_(-zei9) e -" /2 ) ^ = (54) 

with K±(p) = Yia=i(P ~ m a )• This is a "quantization" of the spectral curve [8] K(p) = z + with the 
SW differential dS = p^f. where z = e lx K + {p) and Q(p) = K + (p)K_(p). Even more sophisticated equations 
are needed in the case of toric 1-point function, which is associated with adjoint matter and where one needs 
a quantum deformation of the Calogero system. In the U(2) case it is just given by the Weierstrass function 
potential V(x) — gp(x), while at higher U(n) one has to work with the equation in separated variables. In 
particular, for U(3) [78] 

iij/" + Ul xb" -i(u 2 + 3g(g - l)p{xj)^' - (u 3 + g(g - l)p(x)ui - ig(g - l)(g - 2)p'(x))^ = (55) 

Ui here are the conserved quantities and g is the Calogero coupling constant. 

More detailed comparison of (1) at e 2 = with [57] is rather straightforward, however, one should take 
into account that (49) is not the conventional form of the SW differential for Nf — 2N C implied by the Lax 
formalism for the XXX magnet [8]. In particular, in (49) it is not quite simple to take the limit to the pure 
gauge theory, where [57] provides a nice description. 

It is instructive to study various other limits of formula (1). In fact, the entire variety of limits arising in 
AGT theory is of interest: that of large 0:1,2,3,4, leading from Nf — 2N C case to the pure gauge theory Nf = 
[26]; that of large intermediate dimension a ~ Aq — 7V 2 , associated with Zamolodchikov's asymptotic [37] of 
conformal blocks [27]; that of the large central charge c, giving rise to the ordinary hypergeometric series [28]; 
and the already mentioned limit of e 2 = [56, 57], where the exact BS periods are solutions to the Baxter 
equations. Unfortunately, all these limits are not so easy to take, and each one deserves a separate discussion. 
It deserves emphasizing that study of the DV phases is non-trivial. In the literature there is something known 
about the prepotential expansion in positive powers of S- variables (this is known as the problem of CIV potential 
[52]). What we need now is rather an expansion in negative powers of differences between the 5- variables. This 
is a tedious calculational problem, still relatively straightforward, as we saw above in section 5. However, taking 
limits in this sophisticated procedure is an additional non-trivial step. A check in a rather simple situation of 
= 1 and Nf = 4 — ► Nf = 3,2 in [30] is a nice illustration of the difficulties one needs to overcome. 
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8. Check operators. Higher genus corrections can be analyzed in a similar way, starting from the full 
loop equation, 

p P {zf + (3e 2 V(z)p(z) + (W'(z) + v\z)) P (z) + {fi - l)e 2 ^ = [(w'(z) + v'(z))p(z)] ^ (56) 

The auxiliary potential v(z) — Y^k=o^k zk is used to produce higher multi-resolvents by application of the 
operators V(z) = Y^kLo P^Fi^T' see tne extended review [47] for details. The projector at the r.h.s. of (56) is 
defined 

J+ \ a Z qa I a Z qa 

A reasonable technique to deal with these kind of quantities is provided by the check operators [49] 



For example, 



R{z) = [W'(z)V(z)} + = rzVV(? ) (58) 

z q a 



where 



\4y 2 2y 2 J y y dz 

4y 2 16y \\z z — 1 z — quv J " \z 2 ' (z — l) 2 ' (z — quv) 2 




-^ = ^7.[lz + — + — ) -2(3 + 7—^2+77—^2 ( 6 °) 



Note that the ^4-periods and thus the moduli a arc usually non-affected by higher genus corrections, but these 
additional corrections can change the situation in the present case. 

9. Zpf e k as a T-function. The e-deformation in (1) is in the direction of ei = —£2, where the central charge 
c = 1, thus (3 = 1 and one deals with an ordinary matrix model. It is well known that inclusion of contributions 
of all genera in the matrix model partition function converts the quasiclassical (Whitham) T-function in the 
spherical approximation into the full KP/Toda T-function [46, 79]. Thus we know the meaning of one of the 
two deformations in (1). The other deformation, in the direction of ei with £2 = is known [57] to convert 
the quasiclassical Bohr-Sommcrfcld periods in the formulation of SW theory into the exact BS periods, thus we 
also know the meaning of the other deformation. From the point of view of integrability theory, this implies 
that the ^-function, of which the exact BS periods are monodromies, is the Baker- Akhiczcr function of some 
full (dispersionful) integrable system (Toda, Calogero, Ruijsenaars or magnetic) obtained by switching from the 
Whitham Lax representation to that in terms of differential operators. At the same time, the other deformation 
implies that new time- variables are introduced. Unfortunately, these two deformations are not orthogonal. The 
whole story is a generalization of two-step deformation of the equation p 2 — u = 0: first into the Shroedingcr 
equation Lip = (d 2 — u)ip = and second into the KdV equation, where u becomes a solution to the additional 
hierarchy of equations J-p- = [L k+1 / 2 ) + . This picture can be similar to the vision expressed in [56]. 

10. Conclusion. To conclude, we presented a complete description of Nekrasov functions with arbitrary 
values of parameters £1, £2 and the AGT associated conformal blocks in terms of the exact BS periods and thus 
in the original terms of SW theory. The main ingredient is the SW differential with two deformation parameters 
£1 and £2, and it is given by the matrix model resolvent (1), analytically continued to the DV phase and for j3, 
which can be different from unity. This conjecture still needs to be thoroughly examined and proved, but there 
remain few doubts that it is true. Time is also coming to transfer knowledge from one of the so unified subjects 
to others. As a simple, but still impressive illustration we explained how the Dotscnko-Fateev integrals are now 
fully matched to the generic Virasoro conformal blocks. There are many more applications of this kind to come 
in the close future. 
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